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We have studied the kinetic roughening of an oil-air interface in a forced imbibition experiment in 
a horizontal Hele-Shaw cell with quenched disorder. Different disorder configurations, characterized 
by their persistence length in the direction of growth, have been explored by varying the average 
interface velocity v and the gap spacing b. Through the analysis of the rms width as a function of 
time, we have measured a growth exponent /3 ~ 0.5 that is almost independent of the experimental 
parameters. The analysis of the roughness exponent a through the power spectrum have shown 
different behaviors at short (ai) and long (02) length scales, separated by a crossover wavenumber 
q c . The values of the measured roughness exponents depend on experimental parameters, but at 
large velocities we obtain qi ~ 1.3 independently of the disorder configuration. The dependence of 
the crossover wavenumber with the experimental parameters has also been investigated, measuring 
q c ~ v 0Ar for the shortest persistence length, in agreement with theoretical predictions. 
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I. INTRODUCTION 

The kinetic roughening of growing surfaces is a prob- 
lem of fundamental interest in noncquilibrium statistical 
physics. The interest arises from theoretical, experimen- 
tal, and numerical evidence of scale invariance and uni- 
versality of the statistical fluctuations of rough interfaces 
in a large variety of systems Jj], ||, [|. One candidate 
system is the roughening of a driven interface separating 
two fluids in a porous medium. This problem has im- 
portant practical applications, and has time and length 
scales easily accessible in the laboratory. It allows differ- 
ent possible realizations, depending on the relative vis- 
cosities and wetting properties of the fluids involved ||] . 

One of these possible realizations which has received 
considerable attention in recent years is imbibition, i.e. 
the situation in which a viscous wetting fluid (typi- 
cally oil or water) displaces a second less- viscous, non- 
wetting fluid (typically air) which initially fills the porous 
medium. The motion of imbibition interfaces can be 
spontaneous, i.e. driven solely by capillary forces, or 
forced externally at either constant applied pressure or 
constant injection rate. 

Although there have been many experimental investi- 
gations of the scaling properties of imbibition interfaces 
in the last years, some results, particularly the quanti- 
tative values of scaling exponents, remain controversial. 
The current situation for the case of spontaneous imbibi- 
tion is reviewed in Ref. ||. The situation for the case of 
forced imbibition is summarized in Section [□]. The lim- 
itations of the experiments, in our opinion, arise firstly 
from the lack of precise knowledge of the properties of 
the disorder introduced by the model porous medium, 
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and secondly from the related difficulty in tuning the rel- 
ative strength of stabilizing to destabilizing forces in the 
flow. In the present work we have attempted to avoid 
these two limitations by using a particular model porous 
medium, consisting on a Hele-Shaw cell with precisely 
designed and controlled random variations in gap spac- 
ing. 

In our setup, an initially planar interface becomes sta- 
tistically rough on a mesoscopic scale, as a result of the 
interplay between (i) the stabilizing effects of the vis- 
cous pressure field in the fluid, and the surface tension 
in the plane of the cell, on long and short length scales 
respectively, and (ii) the destabilizing effect of local fluc- 
tuations in capillary pressure, arising from the random 
fluctuations in gap spacing, on short length scales. Al- 
though they have a different physical origin, the role of 
local fluctuations in capillary pressure, in our setup, is 
very similar to the role of wettability defects in the im- 
perfect Hele-Shaw cell introduced by de Gennes ^ , and 
studied by Paterson and coworkers |?J . 

In this paper we present a systematic experimental 
study of forced imbibition of our model porous medium, 
by a wetting silicone oil driven at constant flow rate. 
Since the competing forces in our system are the same as 
in a real porous medium, the physics governing roughen- 
ing is very similar in the two cases. The relative impor- 
tance of viscous forces can be finely tuned by changing 
the injection rate of the invading fluid, and the strength 
of capillary fluctuations can be tuned by adjusting the 
distance between the two glass plates. The sequence of 
photographs in Fig. [l] provides examples of the resulting 
interfaces. 

The outline of the paper is as follows. Section || 
reviews the scaling properties of rough interfaces, and 
their experimental characterization in two-dimensional 
forced imbibition. Section III describes the experimen- 



tal setup, Section IV introduces several parameters, such 
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FIG. 1: Eight consecutive close-up views of the oil-air inter- 
face advancing over the disorder. In each picture the oil is 
driven from bottom to top. The experimental parameters are 
gap spacing b — 0.36 mm, average interface velocity v = 0.08 
mm/s, and disorder SQ 0.40. The time interval between im- 
ages is 0.8 s. 



as permeability and modified capillary number, useful to 
characterize the experiments, and Section ^ explains the 
methodology used in data analysis. The experimental re- 
sults are described i n Se ction VI, and are analyzed and 
discussed in Section VII. The final conclusions are given 



in Section VIII 



II. SCALING OF ROUGH INTERFACES 

To be specific, let us consider a two-dimensional sys- 
tem with cartesian coordinates (x, y), of lateral size L in 
the x direction and of infinite extension in the y direc- 
tion. The interface is driven in the direction of positive y, 
and its position at time t is parametrized by the function 
y = h(x, t). We assume that the interface is initially pla- 
nar, h(x, 0) = 0. Let w(l,t) (the rms interfacial width) 
denote the typical amplitude of transverse excursions on 
a scale I, parallel to the interface, at time t. For the 
complete interface: 



w(L,t)= ((h(x,t) - h(t)) 2 ) 3 



1/2 



(1) 



where h(t) = {h(x,t)) x . The notation (. ..) x represents 
an spatial average in the interval [0, L], where L is the 
system size. The standard picture of kinetic roughening 
is summarized in the dynamical scaling assumption of 
Family and Vicsek S : 



j(L,t)=L a g(t/L z ) 



(2) 



where a is the roughening (or static) exponent and z is 
the dynamic exponent. The scaling function g is: 
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This picture assumes that the lateral correlation length 
of the interface fluctuations increases in time as t x l z , in 
the x direction. For all length scales I within this cor- 
relation length the interface is rough, w(l) ~ l a . As the 
lateral correlation length increases, the interfacial width 
increases correspondingly with time as w ~ t a ' z , which 
defines a growth exponent (3 = a/z. Finally, when the 
lateral correlation length exceeds the system size, L, at 
a crossover time t x ~ L z , the interface fluctuations sat- 
urate. 

One way of measuring a is through the analysis of the 
interfacial width as a function of different system sizes 
L, i.e. w(L) ~ L a , at saturation. From an experimen- 
tal point of view, however, it is not practical to perform 
experiments at different L. Usually it is prefered to mea- 
sure a by analyzing w for different window sizes I, with 
/ ^ L. The scaling of w(l) gives a local roughness ex- 
ponent a/ oc , which is usually identified with the global 
exponent a. This kind of analysis must be performed 
with some caution, however, because it can only pro- 
vide ai oc ^ 1 and hence fails for super-rough interfaces 
(a > 1). Moreover, a number of problems of kinetic 
roughening present intrinsic anomalous scaling, charac- 
terized by different values of local and global exponents 

ii . 

These difficulties can be overcome by analyzing the 
power spectrum of the interfacial fluctuations, defined as 



S(q,t) = (H(q,t)H(-q,t)), 



(4) 
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where 



H(q,t) = [ h (x,t) - h{t) 



(5) 



The notation (...) indicates average over disorder config- 
urations. The mean width w is related to S(q, t) through 



A 



(0) 



where A is the sampling interval in the x direction. For a 
2-d system, the equivalent of the Family- Vicsek scaling 
assumption for the power spectrum reads fl9|: 



S(q,t) = q-^s(qt 1 /> 



where the scaling function is given by 

< 2q+1 for it < 1, 



s(u) 



u 

constant 



for 



u > 1. 



(J) 



(8) 



The power spectrum can give values a > 1 and hence 
is applicable to super-rough interfaces. On the other 
hand, the presence of intrinsic anomalous scaling can be 
detected by a systematic shift of the power spectra com- 
puted at successive time intervals Q. 

The scaling concept has allowed a classification of ki- 
netic roughening problems in universality classes charac- 
terized by different families of scaling exponents. The 
first class is described by the thermal Kardar-Parisi- 
Zhang (KPZ) equation 111]], which provides a local de- 
scription of interfacial roughening in the presence of an 
additive white noise. The KPZ scaling exponents are 
(3 = 1/3 and a = 1/2 in two dimensions. If the non- 
linear term of the KPZ equation is supressed, the result- 
ing equation is known as the thermal Edwards- Wilkinson 
(EW) equation Q, which gives /3 = 1/4 and a = 1/2 in 
two dimensions. When, instead of being purely thermal, 
the noise in the KPZ equation is supposed to depend 
on the interface height h, the equation displays a depin- 
ning transition. At the pinning threshold the interface 
behaviour of the complete equation with the nonlinear 
term ("quenched KPZ") can be mapped to the directed 
percolation depinning (DPD) model, whose scaling expo- 
nents in two dimensions are (3 = a ~ 0.633. Suppressing 
the nonlinearity yields the "quenched EW" equation, for 
which [3 ~ 0.88 and a ~ 1 in two dimensions JL3"[ . 

In the case of imbibition, there are two important is- 
sues which make the above classification of limited ap- 
plicability. The first one is the quenched nature of the 
disorder. It has been argued that for very large driving 
the disorder may still be considered as a fluctuating in 
time (thermal) noise, but in general the disorder must 
be treated as a static, quenched noise. The second issue 
is the nonlocal character of the dynamics, due to fluid 
transport in the cell |l4|, [l5|]. This issue has received 
important attention recently, with the introduction of 
imbibition models which take fluid transport explicitly 



into accoun t, g iving rise to nonlocal interfacial equations 

H 0, H pj . 

The models presented in |l(| [I?], |l8|, |l9| are consistent 
with the well known macroscopic equations of the prob- 
lem (Darcy's law and interfacial boundary conditions). 
They differ in the way the noise is included in the equa- 
tions and in the noise properties. The starting point of 
Ganesan and Brenner |lq| is a random field Ising model. 
The permeability is taken spatially uniform and the noise 
exhibits long range spatial correlations. In the case of 
forced imbibition, and based on a Flory-type scaling, 
the model predicts that the roughness exponent a de- 
pends on the capillary number Ca, with asymptotic val- 
ues a — 3/4 for the smallest drivings and a = 1/2 for the 
largest drivings. The approaches of Dube et al. |l^, 
and Hernandez-Machado et al. [[l7j are both based in 
a conserved Ginzburg-Landau model, where the noise is 
introduced in a fluctuating chemical potential or in the 
mobility, respectively, without long range correlations. 
Dube et al. |H| |l^], however, consider only the case 
of spontaneous imbibition, which could give results very 
different from the case of forced imbibition, specially for 
the growth exponent (3. By numerical integration they 
obtain a = 1.25 and f3 — 0.3. Hernandez-Machado et 
al. [[lTj study the case of forced imbibition, and predict a 
different scaling of the short and long length scales, with 
exponents /3i = 5/6, ati = 5/2 in the former regime, and 
(3 2 = 1/2, OL 2 = 1/2 in the latter. 

A common feature of these models, pointed out in Refs. 
]j~7| , pi |i~9| , is the presence of a new lateral length scale, 
£ c , related to the interplay of interfacial tension and liq- 
uid conservation. For q 3> l/£ c the dominant stabilizing 
contribution is the interfacial tension in the plane of the 
cell, while for q -C l/£ c it is the fluid flow. The interplay 
leads to a dependence of q c (= l/£ c ) on v of the form 



,1/2 



As mentioned in the introduction, the experimental 
characterization of the scaling properties of interfaces 
in two-dimensional forced imbibition remains conflicting. 
Most experiments have been conducted in model porous 
media consisting of air-filled Hele-Shaw cells packed with 
glass beads. In a first experiment of this sort on water-air 
interfaces, carried out by Rubio et al. p0| , a roughen- 
ing exponent a — 0.73 ±0.03 was measured, independent 
of Ca and bead size, for Ca values in the range 10~ 3 - 
10~ 2 . A controversial reanalysis of their data by Horvath 
et al. gave a = 0.91 ± 0.08 |2l|, |2^], which these authors 
compared to a = 0.88 ± 0.08 obtained in their own repli- 
cation of the experiment. In a subsequent work [|23| , 
using glycerol instead of water, Horvath et al. reported 
a clear power law growth regime with a growth exponent 
j3 ~ 0.65, and different values a ~ 0.81 and a ~ 0.49 at 
short and long wavenumbers, respectively, at saturation. 
The last set of experiments of this kind is due to He et al. 
|]i~5| , who explored a very large range of Ca (from 10~ 5 to 
10~ 2 ) and found large fluctuations of the roughness ex- 
ponent in the saturation regime, between 0.65 and 0.91. 
In these experiments w was shown to fluctuate wildly 



4 



during growth, and it was impossible to give a value of 
the growth exponent (3. 

Other experiments have been directed to characterize 
the statistical properties of the avalanches displayed by 
imbibition fronts at sufficiently small Ca. The results 
are also conflicting. The first study EJl was done on 
the same air-glycerol interfaces of Ref. |23fl, and gave a 
power law distribution of avalanche sizes. The second 
one |25f| , on air-water interfaces, reported an exponential 
distribution. 

In our experiments, the possibility of tuning the differ- 
ent competing forces has allowed an accurate measure- 
ment of the growth exponent /3, by enlarging the growth 
regime before saturation. Fine tuning of the forces has 
also allowed to measure the crossover length £ c , and its 
dependence on velocity v. 



III. EXPERIMENTAL SETUP 

In our experiments a silicone oil (Rhodorsil 47 V) dis- 
places air in a horizontal Hele-Shaw cell, 190 x 550 mm 2 , 
made of two glass plates 20 mm thick. The oil has kine- 
matic viscosity v = 50 mm 2 /s, density p = 998 kg/m 3 , 
and surface tension oil-air a = 20.7 mN/m at room tem- 
perature. Fluctuations in the gap thickness are provided 
by a fiber-glass substrate, fixed on the bottom glass plate, 
containing a large number of copper islands which ran- 
domly occupy the sites of a square grid. The height of 
the islands is d — 0.06 ± 0.01 mm. The gap spacing b, de- 
fined as the separation between the substrate and the top 
plate, is set by placing several calibrated spacers on the 
perimeter of the substrate, over the disorder, as shown 
in Fig. ^. We have used gap thicknesses in the range 
6 = 0.16-0.75 (±0.05) mm. 

The disorder pattern is designed by computer and 
manufactured using printed circuit technology. In or- 
der to clearly identify the oil-air interface when it moves 
over the copper islands, the plates have been chemically 
treated to accelerate copper oxidation. Otherwise the 
copper islands are too bright to recognize the interface 
contour. The silicone oil wets perfectly both fiber-glass 
and copper, and no differences in the wetting properties 
due to the oxidation protocol have been observed. 

For each experiment, the glass plates are cleaned with 
soap and water and rinsed with distilled water and ace- 
tone. The fiber-glass plate is cleaned using blotting- 
paper only, which leaves a thin layer of oil to avoid further 
oxidation of the copper surface, ensures a complete wet- 
ting in both fiber-glass and copper, and improves the in- 
terface contrast in the captured images. Next, the fiber- 
glass plate is fixed over the bottom glass plate using a 
thin layer of removable glue. The cell is completed by 
placing the spacers, the upper glass plate, and the o-ring 
around the fiber-glass plate (Fig. g). Finally, the two 
glass plates are firmly clamped together, to ensure the 
homogeneity of the gap spacing and the tightness of cell 
sides closed by the o-ring. 




Fiber-glass 
- plate ( c ) 




SQ SQ-n 



FIG. 2: Sketch of the experimental setup, a) Side view, b) 
Top view, c) SEM image of the copper islands on the fiber- 
glass plate, d) and e) Partial views of the disorder pattern for 
d) squares (SQ, SQ-n) and e) tracks (T). The copper islands 
are the black regions. The blow-ups in d) show the detailed 
structure of the disorder pattern for SQ and SQ-n. 



The oil is injected at constant flow-rate using a syringe 
pump Perfusor ED-2. The syringe pump can be pro- 
grammed to give flow-rates Q in the range 1-299 ml/h 
with less than 2 % fluctuations around the nominal value. 
The oil enters the cell through two wide holes, drilled on 
the top plate near one end of the cell. The other end of 
the cell is left open. To start the experiment with as flat 
an interface as possible, the oil is first slowly injected on 
a transverse copper track, on the fiber-glass plate, which 
is 2 mm ahead the disorder pattern. Next, the syringe 
pump is set to its maximum injection rate until the whole 
interface has reached the disorder (about 3 s later). The 
pump is then set to the nominal injection rate of the ex- 
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FIG. 3: Interface profiles for two types of disorder configura- 
tions, SQ 1.50 (a) and T 1.50 (b), compared with the origi- 
nal digital images. The copper islands are the brighter spots 
of the image background. The experimental parameters are 
b = 0.36 mm and v = 0.08 mm/s in the two cases. 



periment, and t = is defined as the time at which the 
average height of the interface (measured on the images) 
reaches the preset nominal velocity. 

The oil-air interface evolution is monitored using two 
JAI CV-M10BX progressive scan CCD cameras (each 
camera acquiring half side of the cell). The 1/2" CCD 
sensor contains 782 (H) x 582 (V) pixels. In our experi- 
ment we have used an exposure time of 1/25 s in order to 
minimize the illumination. Each camera is equipped with 
a motorized zoom lens Computar M10Z1118MP with a 
focal length in the range 11-110 mm (1:10 zoom ratio). 
The cameras are connected to an Imaging Technology 
PCVision frame grabber installed in a personal com- 
puter. A Visual Basic application controls both cameras 
and stores the images for further analysis. 

The images are taken with a spatial resolution of 0.37 
mm per pixel, a size of 768 x 574 pixels, and 256 grey 
scale levels per pixel. The acquisition is logarithmic in 
time, with temporal increments that vary from 0.33 s to 
180 s. Between 100 and 300 images are taken per exper- 
iment. Because the background is the same for all the 
images captured with the same camera, the interface is 
enhanced by subtracting the first image to all other im- 
ages, and thresholding the result to get a black and white 
contour of the interface. The contour is resolved with 1— 
pixel accuracy using edge detection methods, and data 
are stored for further analysis. This method is automati- 
cally performed, with an error in the interface recognition 
comparable to the width of the oil-air meniscus, about 
one half the gap width. Fig. || presents two examples of 
the digitized interfaces, compared with the original im- 
ages. 



IV. CHARACTERIZATION OF THE 
EXPERIMENTS 

A. Disorder properties 

Three kinds of disorder patterns have been used. Two 
of them are obtained by random selection of the sites of 
a square lattice (Fig. ||(d)). In the first of them (SQ) we 
allow nearest neighbour connections only, leaving next- 
nearest neighbours separated by 0.08 mm. In the sec- 
ond one (SQ-n) both nearest neighbour and next-nearest 
neighbour connections are allowed. The third kind of dis- 
order pattern (T) is formed by parallel tracks, continuous 
in the y direction and randomly distributed along x (Fig. 
§(e)). The filling fraction / (fraction of lattice sites oc- 
cupied by copper) is 35 % in the three cases. However, 
I, the persistence length of the disorder in the direction 
of growth, is very different in the three different disorder 
patterns. This length is measured in the following way: 
(i) We consider every site x of the lattice along the lat- 
eral direction, and measure the average length l{x) of the 
island formed by the connected copper sites at x — Ax, x, 
and x + Ax, where Ax is the lattice spacing, (ii) We av- 
erage l(x) over the lateral direction in the interval (0, L). 
As can be seen in the inset of Fig. |^, I increases by a 
factor 2.4 when changing from SQ to SQ-n, and up to 
the total length of the cell for T. Another characteristic 
of the disorder is the cluster size, which gives the num- 
ber of disorder unit cells of the copper aggregations. The 
statistical distribution of copper clusters is shown in the 
main plot of Fig. [|. 

We have used two values of the lattice spacing in the 
lateral direction, Ax = 1.50 ± 0.04 mm and 0.40 ± 0.04 
mm. From here on we refer to the disorder used in a 
given experiment by SQ, SQ-n, or T, followed by the 
lateral size in mm of the disorder unit cell, 1.50 or 0.40. 



B. Permeability 

To characterize the fluid flow through our Hele-Shaw 
cell with disorder, we have determined the permeability 
of the cell for the different disorder patterns as a function 
of the gap spacing. The experimental set up used for this 
purpose is similar to the one shown in Fig. |2| but the 
injection system has been replaced by a constant pressure 
device. It consists of an oil column of adjustable constant 
height in the range from 200 ± 2 mm to 1000 ± 5 mm. 
The permeability is determined by measuring the oil-air 
interface average velocity for different applied pressures 
(heights of the oil column) and using Darcy's law, 

v = --Vp, (9) 
M 

where v is the interface velocity, k the permeability, fx 
the dynamic viscosity, and Vp the pressure gradient. 
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FIG. 4: Statistical distribution of the copper clusters for a fill- 
ing fraction / = 0.35. Clusters with only nearest neighbour 
contacts (SQ) are represented in dark grey. Clusters with 
nearest and next-nearest neighbour contacts (SQ-n) are rep- 
resented in light grey. The inset shows the persistence length 
of the disorder I as a function of the filling fraction /. The 
solid and dashed lines correspond to the SQ and SQ-n cases 
respectively. The vertical line indicates the filling fraction 
used in our experiments. 

Fig. H shows the results obtained. At large gap spac- 
ings, d/b — ► 0, the disorder has no effect on the fluid flow, 
and the permeability tends to the expected value for an 
ordinary Hele-Shaw cell, ko — b 2 /12, independently of 
the disorder configuration. At very small gaps, d/b — > 1, 
the permeability decreases, and tends to a non-zero value 
k% for d/b = 1, that clearly depends on the disorder con- 
figuration. We have found it convenient to write k\ in 
the form 

ki = d 2 A/12, (10) 

where A is a function that depends on the porosity </> and 
the geometry of the disorder. The simplest functional 
form that interpolates between these two limits can be 
written as: 

The coefficient A can be obtained in general from a fit of 
the permeability to the experimental data. 

In the particular case of disorder T and d/b = 1, an 
analytic expression of A can be derived by recognizing 
that the cell in this case is formed by a parallel array of 
rectangular capillaries. Following Avellaneda et al. [ p6[ , 
in this geometry A is directly the porosity <f>, which in the 
limit d/b = 1 is given simply by <f> = 1 — /. Hence, since 
we have / = 0.35, we obtain A = 0.65. This result fits 
well the experimental data not only in the limit d/b = 1 




FIG. 5: Measurements of the permeability, k, as a function 
of the disorder strength, d/b, for d — 0.06 mm. The dotted 
fine is the permeability of an standard Hele-Shaw cell (with- 
out disorder), and the dashed and solid lines are fits to the 
experimental data for SQ and T disorder configurations, re- 
spectively. The arrows point to the values of d/b for the gap 
spacings used in the experiments. 



but also in the whole range of d/b (solid line in Fig. ||). 
Notice also that the results presented in Fig. [| show that 
there are no important differences between T 1.50 mm 
and T 0.40 mm. This observation generalizes the theo- 
retical result that in the limit d/b — 1 the width of the 
rectangular capillaries does not modify the permeability 
p6[ . For SQ and SQ-n disorder, due to the difficulty of 
finding a general expression for k\ (see Refs. [28]), we 
have fitted k to our experimental results (dashed line in 
Fig. ||) and have obtained a numerical value A = 0.20. 



C. Capillary pressure 

Assuming local thermodynamic equilibrium, the cap- 
illary pressure jump at the interface is given by 

P + ~P~ = c(k|| + k±), (12) 

where k\\ and k± are the main curvatures of the interface, 
in the plane of the cell and perpendicular to it, respec- 
tively. 

Kn varies from for d/b <C 1, to 2/r for d/b ~ 1, 
where r ~ 1.75 mm is the average diameter of the copper 
obstacles. In the range of gap spacings explored we have 
measured curvatures in the range 3.3 x 10~ 3 ^ km ^ 1.15 
mm -1 . Notice that km <C k± in all the range of gap 
spacings used in the experiments. 

In ordinary Hele-Shaw flows (without disorder) k±_ is 
roughly the same in all points of the interface and there- 
fore adds only a constant contribution to the pressure 
jump. In our case, however, when the interface is over 
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the copper islands (gap thickness b — d) we have: 



1 



K±= b VI 



d/b 



(13) 



and when it is over the fiber-glass substrate (gap thick- 
ness b): 



(14) 



assuming complete wetting. 

This difference in curvature, given by 



d/b 
b \l-d/b 



(15) 



makes the interface experience a capillary instability 
when it passes from one gap spacing to the other. In 
the range of gap spacings studied, 0.16 < b ^ 0.75 mm, 
we get 7.5 ^ ^ ^ 0.23 mm -1 . We have verified that for 
gap spacings b > 2 mm, which correspond to ^ < 0.031 
mm , the fluctuations in capillary pressure are no more 
sufficient to roughen the interface appreciably. 



D. Modified capillary number 

Once the permeability of the cell has been character- 
ized, we can introduce a dimensionless number to de- 
scribe the relative strength of viscous to capillary forces. 
The simplest number that relates viscous and capillary 
forces is the capillary number Ca= fj,v/a. In order to 
account for the properties of the disorder, which are not 
contained in the previous definition of Ca, it is custom- 
ary to introduce a modified capillary number. For an 
ordinary Hele-Shaw cell (without disorder), the modified 
capillary number, which we call Ca*, comes out from the 
dimensionless form of the Hele-Shaw equations p9[ : 



(V = ( a • 12 | | 



(16) 



To define a modified capillary number Ca' for our par- 
ticular cell with disorder, we consider on one side the 
average viscous pressure drop across the cell, given by 



Sp v 



L ■ \V Pv 



(17) 



with \7pvis = —vji/k (Darcy's law). L, the cell width, 
provides the macroscopic length scale. On the other side, 
a measure of the capillary pressure drop is given by 



Spcap = o- 



b-d 



(18) 



The last contribution accounts for the curvature of the 
interface in the plane of the cell, and is relevant only 
when d/b — > 0, i.e. when the destabilizing role of the 
disorder vanishes. 
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FIG. 6: Ratio of the modified capillary number for the cell 
with disorder, Ca', to the modified capillary number for an 
ordinary Hele-Shaw cell, Ca*, as a function of the disorder 
strength, d/b. The dashed and solid lines correspond to dis- 
orders SQ 1.50 and T 1.50. The dotted line represents the 
asymptotic limit Ca'=Ca* . The arrows point to the values of 
d/b for the gap spacings used in the experiments. 



Defining Ca' = 5p vis /Sp cap and using (|ll|), (|l7|) and 
@, we get 



Ca' = Ca • 



12L 



*(>-<i-a">!) (m + x 

(19) 

The ratio Ca'/Ca* as a function of d/b is shown in Fig. 
^ for disorders SQ and T. When d/b — > the destabiliz- 
ing role of the disorder is negligible, as expected, and 
Ca'/Ca* tends to 1 (ordinary Hele-Shaw cell). As d/b in- 
creases, the increasing strength of the disorder is manifest 
in the progressive decrease of Ca'/Ca*. It is important to 
notice that our definition of Ca' is not valid for d/b = 1, 
because the flow would be essentially different when the 
free gap disappears. In the range of d/b explored in our 
experiments the permeability of the cell remains always 
close to that of an ordinary Hele-Shaw cell (Fig. ||), and 
the decrease of Ca'/Ca* is essentially due to the capillary 
forces associated with the menisci in the z direction. 

In the range of gap spacings experimentally explored, 
0.16 ^ b ^ 0.75 mm, the ratio between the free gap b — d 
and the total gap b varies between 62% and 92%. These 
large ratios, combined with the high viscosity and small 
surface tension of the silicone oil compared with other 
fluids (i.e. water), are responsible for the large values 
of Ca' in our experiments, in comparison to the values 
reported for pure porous media. Considering for example 
the disorder T 1.50 mm, Ca' varies from a value 1.33 
(for b — 0.16 mm and the minimum interface velocity, 
0.04 mm/s) to a value 17.0 (for b = 0.75 mm and the 
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maximum interface velocity, 0.4 mm/s). 



V. DATA ANALYSIS 

The wetting of the lateral gap spacers by the invad- 
ing oil changes the physics at the two ends of the in- 
terface. To minimize this disturbance, which is partic- 
ularly important at large gap spacings, we have disre- 
garded 8 mm at each side of the cell, thereby reducing 
the measured interface from 190 to 174 mm in the x 
direction. The final number of pixels of the interface 
after this correction is reduced from 515 to 470. For 
data analysis convenience these pixels are converted into 
N = 512 equispaced points, through a linear interpola- 
tion. The final spacing between two consecutive points 
is A = L/N = 0.34 mm . Although the linear interpo- 
lation introduces an artificial resolution 7 % larger than 
the resolution of the original image, the increase does not 
affect the final analysis. 

The interfaces measured at the smallest gap spacings 
and velocities may present overhangs. These multivalua- 
tions, rather exceptional, have been eliminated by taking 
for each value of x the corresponding largest value of 
h(x). 

In addition, we have forced periodic boundary condi- 
tions for h(x,t) by subtracting the straight line connect- 
ing the two ends of the interface. This procedure is well 
documented in the literature of kinetic roughening [j30| . 
The linear correction imposed to the interfaces has signif- 
icant effects on the power spectrum. The Fourier spec- 
trum of an interface which is discontinuous at the two 
end points is dominated by an overall behaviour of the 
form q~ 2 . Forcing periodic boundary conditions elimi- 
nates the overall slope -2 |3l]]. Moreover, the analysis of 
w(l) is insensitive to the linear correction of the interface 
(except for values of I comparable to the system size) and 
gives results consistent with those obtained from S(q, t) 
only for interfaces with periodic boundary conditions. 

Since the maximum width of the meniscus is only one 
half the gap width (in conditions of complete wetting) the 
interface can be considered one-dimensional at the length 
scale of the copper islands. The resolution of the interface 
fluctuations in the y direction is ±1 pixel at a given point. 
Nevertheless, the measurement of the global interfacial 
width w is much more precise because the width is an 
average over the N = 512 points of the interface. 

Given that, according to our choice of the time ori- 
gin, the whole interface is already inside the disorder at 
t = 0, w(0) ^ 0. For this reason, we have decided to 
characterize the interface fluctuations by the subtracted 
width W, defined by W 2 {t) = w 2 {t) - w 2 (0) g, |§. No- 
tice that, since the data analysis is based on power law 
dependencies, short times are very sensitive to the defi- 
nition of t = and to the value w(0). After analyzing 
the data for different definitions of t = 0, and checking 
the influence of subtracting u>(0), we have found that 
the analysis based on the subtracted width is the most 



objective and less sensitive to the details of the exper- 
imental procedure. The error bars shown in the W(t) 
plots indicate the dispersion of the different individual 
experiments with respect to the average curve plus the 
uncertainty in the determination of t = 0. 

Finally, the crossover time t x has been measured on 
the W(t) log-log plots as the time when the power law 
with slope j3 crosses the horizontal straight line that cor- 
responds to the average value of the interfacial width at 
saturation, W s . 

VI. EXPERIMENTAL RESULTS 

The parameters explored in our experiments are sum- 
marized in Table |. The minimum velocity selected in 
the experiments, v — 0.04 mm/s (which will be taken as 
reference unit for interfacial velocities), has been chosen 
to ensure that the interface is always single-valued for 
gap spacings b > 0.36 mm. We have used three differ- 
ent disorder realizations for SQ and SQ-n and four for T. 
For each disorder realization and each set of experimen- 
tal parameters (v, b) we have carried out three runs for 
SQ and SQ-n, and two runs for T. 

A. Variable velocity 

1. SQ 1.50 mm 

We have started the study exploring five different ve- 
locities, from V to 10 V, using an intermediate gap spac- 
ing b = 0.36 mm. We have not selected velocities lower 
than V to avoid overhangs and trapped air. An example 
of the interfaces obtained at different velocities is pre- 
sented in Fig. 0. All the experiments shown correspond 
to the same disorder configuration, and the average po- 
sition of the selected interfaces is also the same for the 
different velocities. The most noticeable aspect of the se- 
quence of interfaces is that the short length scales are not 
affected significantly by the velocity, as opposite to the 
long length scales, which become progressively smoother 
as the velocity increases. 

Fig. H shows the W(t) plot for four different velocities. 
We have omitted the curve for 7V because it is too close 
to the 10 V curve. All curves display a power law growth 
regime with approximately the same growth exponent 
(3 = 0.47 ± 0.04, independent of the velocity. Saturation 
times and saturation widths depend clearly on the veloc- 
ity, both decreasing at increasing velocities. This general 
behaviour is not so clear for the two smallest velocities, 
for which the corresponding t x and W s are practically the 
same. This can be a consequence of the fact that the in- 
terface at very low velocity gets locally pinned at the end 
of the copper obstacles, hindering an arbitrary deforma- 
tion of the interface. The strong fluctuations of W that 
can be observed both during growth and at saturation 
are also remarkable. These fluctuations have two origins: 
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Disorder 


Gap spacing (mm) 


interlace velocity (V = 0.04 mm/sj 


SQ 1.50 


0.16,0.23,0.36,0.56,0.75 


2V 


SQ 1.50 


0.36 


V, 21/, 51/, 71/, 10 V 


SQ-n 1.50 


0.36 


V, 21/, 10 V 


T 1.50 


0.16,0.36,0.56,0.75 


21/ 


T 1.50 


0.36 


1/, 2V, 51/, 10 V 


SQ 0.4 


0.36 


21/ 


T 0.4 


0.36 


21/ 



TABLE I: Summary of the parameters explored for each kind of disorder. 
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FIG. 7: Sequence of interfaces at different velocities, for b — 
0.36 mm and disorder SQ 1.50, corresponding to Ca'= 1.85, 
3.69, 9.24, and 18.5. The disorder realization is the same in 
the four sequences. The reference velocity is V — 0.04 mm/s. 
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FIG. 8: Interfacial width, W[t), for disorder SQ 1.50, gap 
spacing b = 0.36 mm, and four different velocities. The 
straight line drawn to guide the eye is the power law W ~ 
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the intrinsic metastabilities in systems with quenched-in 
disorder fl5|| , and to a lesser degree the inhomogeneities 
in the gap thickness caused by small deformations in the 
fiber-glass substrate or small differences in the height of 
the copper obstacles. To smooth out the fluctuations it 
would have been necessary to increase the number of dis- 
order configurations and runs to a number that would be 
prohibitive. Nevertheless, the largest fluctuations usually 
appear deep in saturation and do not change the results 
presented here. 

Fig. ^| shows the analysis of the interface fluctuations 
through the power spectrum. An example of the tem- 
poral evolution of the spectrum, which corresponds to 
the experiments at velocity 2V, is presented in Fig. ||(a). 
Since larger length scales become progressively saturated 
with time, the spectrum at short times displays a power 
law decay for large q only. As time increases the power 
law extends to smaller q until a second power law be- 
haviour emerges, which reaches the smallest q, corre- 
sponding to the system size L, at saturation. The first 
power law regime (large q) will be characterized with a 
roughness exponent ot\, and the second one (small q) 
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with a roughness exponent a<i- In the range of param- 
eters explored we have always observed a.\ > a?,. It is 
interesting to notice that the first regime shows a slightly 
time-dependent behaviour. a% is progressively smaller 
as time increases, and reaches a constant value at satu- 
ration. This temporal dependence disappears at cither 
high interface velocities or large gap spacings. 

The power spectrum at saturation for five different ve- 
locities is shown in Fig. ^(b) . The different behaviour at 
short and long length scales is visible in all cases, but the 
crossover from one regime to another (characterized by 
the wavenumber q c at the crossing point of the two power 
laws) increases systematically with increasing velocities. 

The exponent ai (large q) increases with the interface 
velocity, from a value a± = 0.9 ±0.1 (slope = —2.8 ±0.2) 
at velocity V to a value a\ — 1.3±0.1 (slope = — 3.6±0.2) 
at velocity 10 V. We have observed that the variation is 
not linear, but tends to the limiting value a\ — 1.3 at 
large velocities. 

The exponent a 2 (small q) decreases with the interface 
velocity, varying from a 2 — 0.6 ±0.1 (slope = —2.2 ±0.2) 
at velocity V to a 2 = 0.1 ± 0.1 (slope = -1.2 ± 0.2) at 
velocity 10 V. This exponent seems to be insensitive to 
the velocity at low velocities. 

Fig. |9|(c) shows the alternative analysis carried out to 
determine the roughness exponents, using W(l). Notice 
that because this is a local analysis it is limited to a 
maximum value a = 1. At short length scales we obtain 
Qi = 0.90 ±0.1 for all the velocities. The value a.\ — 1 
that we should have obtained for v > 2V is unreach- 
able because the fluctuations in W(l) contribute always 
in the direction of decreasing a%. At long length scales 
the values obtained for a 2 are consistent with the values 
measured from the power spectrum. 



2. Effect of increasing the persistence of the disorder 

We have carried out a series of experiments oriented to 
explore how different kinds of disorder patterns affect the 
interfacial dynamics, and the possible universality of our 
results. The gap spacing and velocity used in all cases 
are b = 0.36 mm and 2V. 

When the correlation length of the disorder in the y di- 
rection is increased (changing the disorder pattern from 
SQ to SQ-n and then to T), we observe important dif- 
ferences in the behaviour of both W(t) and S(q,t). The 
correlation length of the disorder in the y direction is 
quantified through the average extent of the disorder cells 



(a) 



in the y direction, I, introduced in Sec. IV A. For a basic 
disorder cell of size 1.50 mm, I ~ 2 mm for SQ, I ~ 5 mm 
for SQ-n, and I — oo for T. The immediate consequence 
of a progressively larger correlation length is the persis- 
tence of the capillary forces for longer times. The effect 
on W{t) is that the large fluctuations observed during 
growth disappear as I increases. In the extreme limit of 
disorder T, W(t) is a neat power law with an exponent 
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FIG. 9: Determination of the roughness exponents for ex- 
periments at five different velocities, with 6 = 0.36 mm and 
disorder SQ 1.50. a) Example of the temporal evolution of the 
power spectrum for velocity 2V (Ca'=3.69). b) Power spec- 
tra at saturation for the five velocities studied. The curves 
have been shifted in the vertical direction for clearness, c) 
Analysis of the local width, W(l). The vertical lines in (a) 
and (b) give the value of q associated with the spatial scale of 
the disorder. The tilted line crossing the power spectra in (b) 
fits the values of the crossover wavenumber q c at the different 
velocities. The other straight lines are power law fits from 
which the roughness exponents can be deduced. 
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= 0.52 ± 0.03, as can be seen in Fig. JL0|(a). Other rele- 
vant aspects are that the saturation times can be clearly 
identified and the fluctuations in saturation have a mi- 
nor amplitude. Whereas the important fluctuations of 
W(t) both during growth and saturation, together with 
substantial differences between different runs and disor- 
der configurations, are characteristic of the experiments 
with SQ, for T the experiments show that different runs 
with the same disorder configuration lead to practically 
identical results, and differ only slightly when the disor- 
der configuration is changed. 

The increasing correlation length in the direction of 
growth has important consequences in the power spec- 
trum. On one side, the observation of a smaller rough- 
ness exponent a\ as time increases for disorder SQ (Fig. 
^(a)), is more evident as we go to SQ-n and T. The expo- 
nent a.\ obtained for different disorders, however, tends 
to a single value at large gap spacings and high velocities 
(weak capillary forces). On the other side, the spectra 
at saturation shown in Fig. [H)](b), clearly indicate that 
a progressively larger I reduces the measured roughness 
exponents, from ol\ = 1.0 ± 0.1 to a± — 0.5 ± 0.1, and 
from ot2 = 0.6 ± 0.1 to ol-i = 0.2 ± 0.1, as the disorder 
changes from SQ to T. This trend is confirmed by exper- 
iments using SQ 0.40 (I ~ 0.6), where the evolution of 
a\ with time is not observed due to the short I and, at 
saturation, gives exponents a\ ~ 1.1 and a2 — 0.5. 

Another interesting modification that can be intro- 
duced in the disorder is the size of the basic disorder cell. 
In addition to our usual experiments with lateral size 
1.50 mm, we have made experiments with disorder T of 
lateral size 0.40 mm. The W(t) curves for both cases are 
shown in Fig. [l0|(c). Both curves give approximately the 
same growth exponent, around f3 = 0.50. The small dif- 
ferences in slope around this value can be attributed to 
the reduced number of disorder configurations studied. 
The most noticeable difference is the saturation width, 
which is smaller for the small disorder size. The power 
spectrum at saturation for the two cases is shown in Fig. 
|l0|(d). Notice that in the case of T 0.40, the measured 
exponents u\ and ui increase considerably, and approach 
the values measured for SQ or SQ-n 1.50. 



3. T 1.50 mm 

Here we present a systematic study of the disorde r 
T 1.50, in the same line as for SQ 1.50 (Sec. |VIAl| ). 
The disorder T 1.50 ensures that the interface is always 
single-valued at any gap spacing and velocity, which al- 
lows exploring a wide range of experimental parameters. 

Fig. [ll] shows an example of the interfaces obtained in 
the experiments with T 1.50 at different velocities and 
fixed gap spacing. As we observed for SQ 1.50, the in- 
terfaces become progressively flatter as the velocity is in- 
creased. However, there are interesting differences. The 
first one is the possibility to explore the regime of very 
low velocities, thanks to the fact that the continuity of 



the copper tracks makes local pinning impossible. How- 
ever, there are two reasons for not selecting velocities 
lower than V . First, we have observed that for v < V/4 
the oil fingers become so elongated that the interface 
pinches off at some distance behind the finger tip. And, 
second, because the oil tends to advance preferentially 
over the copper tracks, at very short times the oil over 
the fiber-glass recedes (due to mass conservation) and 
parts of the interface reach again the transverse copper 
tracks at the beginning of the cell. This limits the veloc- 
ity to a minimum value V for this gap spacing. 

As discussed before, having tracks impedes that the 
interface gets locally pinned. As a result, the metastabil- 
ities observed for disorder SQ are not present here. Most 
of the fluctuations of W(t) disappear, as shown in Fig. 
|l2| for the four velocities studied. All the curves have ap- 
proximately the same growth exponent (3 = 0.52 ± 0.03. 
The transition from the growth regime to the saturation 
regime is rather smooth and makes it difficult to measure 
the values of t x and W s . Although two identical runs 
(with the same disorder configuration) give almost iden- 
tical W{t) curves, different disorder configurations give 
curves with slightly different (3, t x , and W s . In the for- 
mer case the values vary between (3 = 0.48 and (3 = 0.55. 
The fluctuations observed at saturation are attributed to 
small inhomogeneities in the gap thickness, which affect 
particularly the experiments at the lowest velocities. Al- 
though the amplitude of the fluctuations decreases when 
we increase the number of disorder configurations, the 
almost perfect reproducibility of experiments with iden- 
tical disorder configuration makes it difficult to have good 
statistics using a reduced number of disorder configura- 
tions. 

The analysis of the power spectrum for the four veloc- 
ities is shown in Fig. [L3|. The deterministic character of 
the interfacial growth makes the interfacial fluctuations 
not self-averaging. This leads to considerable fluctua- 
tions in the power spectrum, that cannot be completely 
reduced without using a large number of disorder con- 
figurations. Qualitatively, the power spectrum presents 
the same behaviour observed for SQ, with two different 
regimes separated by a crossover point q c . The rough- 
ness exponent ot\ increases with the velocity, varying 
from ai = 0.5 ± 0.1 (slope -2.1 ± 0.1) for velocity V 
to at = 1.2 ±0.1 (slope -3.4 ±0.1) for velocity 10 V. We 
have observed that at high velocities a± reaches the same 
limiting value as it was measured for SQ 1.50. 



B. Variable gap spacing 

1. SQ 1.50 mm 

Variations in the gap thickness modify the strength of 
the capillary forces in the vertical direction, which are 
destabilizing, in relation to the viscous forces and to the 
capillary forces in the plane of the cell, which are stabi- 
lizing. The sequence of interfaces at different b and equal 
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FIG. 10: Interfacial width W(t) (a) and power spectra S(q,t) (b) for three kinds of disorder with unit cell of size 1.50 mm. c) 
and d) show the same magnitudes for disorder T of two different sizes, 0.40 and 1.50 mm. The spectra have been shifted in the 
vertical direction for clearness. The vertical lines give the value of q associated with the spatial scale of the disorder. In all the 
experiments b = 0.36 mm and v = 2V (Ca'= 3.25 for T, and 3.69 for SQ and SQ-n). 



velocity (2V) of Fig. [l4|, indicate that roughening at large 
scales is inhibited for the smallest gap spacing, b = 0.16 
mm, and the interface remains globally flat. Notice that 
the oil-air interface adapts perfectly to the disorder con- 
figuration. The gap spacings must be increased to values 
larger than b — 0.23 mm to observe all scales contributing 
to the roughening process. 

The results of W(t) for different gap spacings are pre- 
sented in Fig. |l5|(a). The results show that the satura- 
tion times and widths are practically independent of gap 
spacing. This is a demonstration that increasing the gap 
thickness is not equivalent to increasing the velocity, al- 
though both changes result in larger Ca'. The growth 
exponent (3 at the smallest gap spacings centers around 
j3 = 0.40 but is subjected to a large error bar due to 



a significant presence of multivaluations in the interface. 
As b increases there is a clear tendency towards a well de- 
fined power law, over three orders of magnitude in time, 
resulting in an exponent (3 — 0.50 ± 0.02. 

Although all the experiments start with the same ini- 
tial condition, the transition from an interface almost 
flat, with W ~ 0, to a set of almost parallel W(t) curves, 
suggest that the growth at very short times is strongly 
dependent on gap thickness. 

The different role of velocity and gap spacing is also 
apparent from the power spectrum of the interfaces at 
saturation, shown in Fig. |l5|(b). For the smallest gap 
spacing, b — 0.16 mm, the spectrum is nearly flat at small 
g, reflecting that the largest scales do not grow (Fig. [u]) . 
The only contribution to this part of the spectrum comes 
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FIG. 11: Sequence of interfaces for disorder T 1.50 and gap 
spacing b = 0.36 mm, at different velocities, corresponding to 
Ca'=1.62, 3.25, 8.13, and 16.3. The disorder pattern is shown 
in grey. 
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FIG. 12: Interfacial width, W(i), for disorder T 1.50, gap 
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from the memory of the global shape of the initial front, 
which, at the smallest gap spacings, cannot be prepared 
perfectly flat. The power law dependence at larger q 
extends from q c ~ 0.15 to q — 0.67 (which corresponds 
to the lateral size of the disorder cell). As b increases, the 
nearly q-independent behaviour at small q changes to a 
power law behaviour (exponent ai), and the crossover 
value q c shifts to smaller q. For gap spacings as large 
as 6 = 0.75 mm this shift leads to a single power law 
dependence, that could be strongly affected by finite- 
size effects. The measured exponents in the short length 
scale regime increase with b, from a.\ = 0.9 ± 0.1 (slope 
-2.8 ± 0.2) to a t = 1.1 ± 0.1 (slope -3.2 ± 0.2). 



FIG. 13: Determination of the roughness exponents from the 
power spectra for experiments at four different velocities, with 
b = 0.36 mm and disorder T 1.50. The curves have been 
shifted in the vertical direction for clearness. The vertical 
line gives the value of q associated with the spatial scale of 
the disorder. The other straight lines are power law fits from 
which the roughness exponents can be deduced. 



2. T 1.50 mm 

The gap spacing has been varied in the range 0.16 ^ 
b ^ 0.75, keeping the velocity fixed at 2V. Fig. [l(] show 
a sequence of the temporal evolution of the interfaces for 
four different gap spacings. The most remarkable aspect 
is the significant difference in the shape of the interfaces 
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FIG. 14: Sequence of interfaces at velocity 2V and disor- 
der SQ 1.50, at different gap spacings corresponding to, from 
smaller to larger gap spacing, Ca'= 3.65, 3.70, 3.66, and 3.62. 
The disorder realization is the same in the four sequences. 
The reference velocity is V — 0.04 mm/s. 
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FIG. 15: Experimental results at velocity 2V and disorder SQ 
1.50 at five different gap spacings. a) Interfacial width W(t). 
A power law with slope 0.5 has been drawn to guide the eye. 
b) Power spectra at saturation (t = 6000 s). The curves 
have been shifted in the vertical direction for clearness. The 
vertical line gives the value of q associated with the spatial 
scale of the disorder. 



between b = 0.16 mm and the other gap spacings. For 
b = 0.16 mm the capillary forces are strong enough to in- 
hibit correlations between neighbouring tracks. Similarly 
to the experiments with SQ 1.50 and the smallest gap, we 
do not expect saturation of scales larger than the average 
lateral size of the tracks. For b ^ 0.36 mm the behaviour 
is different: the interfacial height is correlated from the 
size of the basic disorder cell up to the system width. 
The interface at long length scales adopts a final shape 
at saturation identical for the three gaps, dominated by 
the spatial distribution of the disorder pattern. At short 
length scales the amplitude of the fingers depends on the 
relative strength of the capillary forces, tuned by the gap 
spacing. 



Fig. [i~7|(a) shows the W{t) curves for the four gap spac- 
ings. The slopes are practically the same in all four cases, 
the saturation times are also very similar, and the satu- 
ration widths are progressively larger as the gap spacing 
is reduced. The fact that the interfacial fluctuations with 
b = 0.16 mm reach saturation is a consequence of driving 
the system with a finite velocity that impedes an infinite 
deformation of the interface. 

In Fig. |l7](b) we represent the spectra for the four ve- 
locities at saturation. The fact that for b = 0.16 mm the 
interfacial height is uncorrelated for distances larger than 
the average lateral size of the tracks (2.3 mm, correspond- 
ing to q = 0.43) is clearly observed, since for q < 0.4 we 
observe a plateau in the power spectrum. The only hint 
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FIG. 16: Sequence of interfaces at velocity 2V and disorder T 
1.50, at different gap spacings corresponding to, from smaller 
to larger gap spacing, Ca'= 2.67, 3.25, 3.38, and 3.41. 
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FIG. 17: Experimental results at velocity 2V and disorder T 
1.50, at four different gap spacings. a) Interfacial width W(t). 
A power law with slope 0.52 has been drawn to guide the eye. 
b) Power spectra at saturation (t = 6000 s). The curves 
have been shifted in the vertical direction for clearness. The 
vertical line gives the value of q associated with the spatial 
scale of the disorder. 




of a power law regime is at large q, from q > 0.4 to the 
limit given by the disorder size, an interval too short to 
measure any roughness exponent. The other gap spac- 
ings show qualitatively the same behaviour observed for 
SQ. The large q saturate with an exponent ai that is 
progressively larger as we increase the gap spacing, vary- 
ing from ax = 0.5 ± 0.1 (slope -2.1 ± 0.2) for b = 0.36 
mm to ai = 0.9 ± 0.1 (slope -2.9 ± 0.2) for b = 0.75 
mm. The regime characterized by ai is only identifiable 
for b = 0.36 mm. For the other gap spacings, our results 
indicate that we have a unique power law over all length 
scales. Here again, however, the results for the largest 
gap spacings could be severely influenced by finite size 
effects. 



VII. ANALYSIS AND DISCUSSION 

The variation of the roughness exponents with velocity 
and gap spacing, for SQ and T disorder configurations, is 
summarized in Fig. We observe that a\ and ai have 
opposite behaviours as the interface velocity increases. 

The exponent a% , characteristic of roughening at short 
length scales, increases with v. It becomes larger than 
1 (super-rough interface) already at moderate veloci- 
ties. To understand this behaviour, we refer to Fig. ^|(a), 
where it is shown that a\ is large at very short times, 
and decreases progressively as time goes on, at velocity 
2V. The initial super-roughness is due to the initially 
local dynamics of the interface fluctuations as the inter- 
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FIG. 18: Variation of the roughness exponents with interface 
velocity (top) and gap spacing (bottom), for SQ 1.50 (circles) 
and T 1.50 (triangles). Solid and open symbols correspond to 
ai and a-i , respectively. The region of super-rough interfaces 
(a > 1) is shown in grey. The lines are drawn to guide the 
eye. 



face gets in contact with the disorder for the first time. 
At sufficiently large velocities, however, the subsequent 
decrease of a\ is not observed, and the initial super- 
roughness gets frozen in. The reason is that the satu- 
ration time t x is comparable to the average time spent 
by the interface to go across the distance I (the average 
length of the disorder in the y direction), estimated as 
t = I /v. For example, for SQ 1.50, b = 0.36 mm and 
v = 10 V, we have t x = 5 s and t = I/u~6s. This ex- 
plains also that a\ takes similar values for SQ and T at 
high velocities, since whenever t x ~ r the continuity of 
the disorder in the y direction is irrelevant. In particular, 
we have evidence that ot\ becomes identical for SQ and T 
at very large velocities. The interface is saturated at all 
time, making the interfacial dynamics in SQ equivalent 
to an average of the dynamics over a number of T config- 
urations. For gap spacings around b — 0.36 mm, we have 
observed that this condition is satisfied for v > 13V. An 
identical value a\ ~ 1.3 has been obtained for both SQ 
and T by performing experiments at 30V" (t x <t). This 
value of a\ is distinctively larger than previous experi- 
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FIG. 19: Saturation time t x as a function of velocity v (top) 
and gap spacing b (bottom), for SQ 1.50 (circles) and T 1.50 
(triangles). The straight lines have been drawn to guide the 
eye. 



mental results, and coincides within error bars with the 
value 1.25 found numerically in Ref. jl9) taking account 
of the nonlinearities. A super-rough interface (a > 1) is 
also obtained by scaling analysis in jlTj , but with a larger 
exponent. 

The exponent cv2 (long length scales) measured in our 
experiments decreases with v. Its actual behaviour, how- 
ever, could be obscured by finite size effects. The reason 
is that nonlocal effects on the interfacial dynamics (due 
to the fluid flow from the reservoir) become progressively 
more significant as v increases. For experimental reasons 
we have not studied finite size effects in any systematic 
way. However, we do have checked that experiments per- 
formed at velocity V and system size L/2 give the same 
exponent ai than those at velocity 2V and system size 
L. The simultaneous dependence of ai on velocity and 
system size can be estimated from W(t) for both SQ 
and T. We have observed (Figs. [n](a) and ^(a)) that 
the saturation time i x and the saturation width W s are 
power laws of the velocity v, of the form t x ~ v~ s and 
W s ~ v" 7 , with 6 ~ 2.5 and 7 ~ 1.1 for SQ, and S ~ 2.3 
and 7 ~ 1.2 for T. Combining these relations with the 
FV scaling assumption (||) , and taking into account that 
the saturation times and widths correspond to the long 
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FIG. 20: Saturation width W s as a function of velocity v (top) 
and gap spacing b (bottom), for SQ 1.50 (circles) and T 1.50 
(triangles). The straight lines have been drawn to guide the 
eye. 



time regime (long length scales), we get: 

. log v 



a 2 



-(36 



logL' 



(20) 



with 7 = (38. Given that the system size L has been 
kept fixed in the experiments, the relation ( |20| ) predicts 
that the roughness exponent ai obtained from the ex- 
periments will depend on v. The validity of this result is 
confirmed in Fig. where the values of a 2 are plotted 
vs. log v. 

From the observation that the difference between ati 
and a 2 decreases as v is reduced (Fig. |l8| ), and taking 
into account that capillary forces become also more im- 
portant, we argue that in the limit of very low velocities 
the power spectrum would display an unique power law 
extending over all scales with a roughness exponent in 
the interval 0.6 — 0.9, compatible with the observations 
in pi |20] M, pi pi. This limit is unreachable for SQ 
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FIG. 21: Roughness exponent ai as a function of logu for SQ 
1.50 (circles) and T 1.50 (triangles). The straight lines have 
been drawn to guide the eye. 



because at velocities v < V/2 the interface gets locally 
pinned and develops overhangs. In the case of T there is 
no pinning or multivaluation, but the interface stretches 
to such a point that easily pinches off at long times, be- 
fore saturation. These and other questions, present only 
in the regime of flows dominated by capillary forces, are 
studied in Ref. |||. 

Concerning the influence of b, the gap spacing, two 
main effects must be taken into account. First, the mo- 
bility is proportional to b 2 and hence increases with 6, 
enhancing the stabilizing effect of the viscous pressure 
field (on long length scales) and the interfacial tension in 
the plane of the cell (on short length scales) . Second, in- 
creasing b weakens the destabilizing effect of the capillary 
forces induced by the disorder, which act on short length 
scales. The overall result is that the curves W(t) are 
shifted towards smaller widths as b increases, as shown 
in Figs. |l|(a) and |l7|(a). The shift is such that W ■ b col- 
lapses the curves for different gap spacings into a single 
one, in both SQ and T cases. The fact that large val- 
ues of b makes the presence of the disorder irrelevant can 
be observed through the values of the roughness expo- 
nents a, which show a tendency to converge into a single 
value a ~ 1 at large gap spacings and for any disorder 
configuration. 

Another observation is that the saturation regime is 
reached always at almost the same saturation time t x , 
independently of the gap spacing, as Fig. |l^(b) reveals. 
The interfacial width at saturation, however, has a dif- 
ferent origin depending on gap spacing. For small b the 
width is due to fluctuations at the shortest length scales, 
while for large b it is due to fluctuations at longer length 
scales. This can be observed on the interfaces and is 
clearly reflected on the power spectrum (Figs. f[5](b) and 
|l7|(b)). The passage from short to long length scales is 
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continuous with b. It can be shown that replacing the 
variable q by the dimcnsionless variable q ■ b, produces 
the collapse of the power law region of the power spectra 
in Fig. ||(b). 

Table || shows the main results obtained for the rough- 
ness exponents a% (short length scales) and a-i (long 
length scales) at different disorder configurations and 
drivings. For small Ca', capillary forces are dominant 
at all scales, and the dynamics is very sensitive to the 
disorder configuration. For SQ and SQ-n, the interfaces 
get locally pinned, and the measured exponents are close 
to those obtained in DPD or in experiments where cap- 
illary forces are dominant Jl5|, [|l], |22|, ^3). In the 
limit of persistent disorder, the nature of the disorder im- 
pedes pinning, but the effect of the destabilizing capillary 
forces combined with the correlations between neighbor- 
ing tracks leads to a new regime that can be described 
using the anomalous scaling ansatz fl33| . This new regime 
also extends to the region of moderate Ca' for T disor- 
der. For moderate Ca', viscous forces are dominant at 
long length scales and two clear regimes separated by 
a crossover wavenumber q c can be characterized. For 
SQ and SQ-n we get roughness exponents ai > 1 and 
«2 — 0.5 — 0.6. For T, we get a qualitatively similar be- 
haviour, but obtaining lower values of the exponents due 
to finite size effects at long length scales and the domi- 
nant capillary forces at short length scales. For large Ca', 
the viscous forces cause that the initial super-roughness 
at short times and short length scales gets frozen in, ob- 
taining the same roughness exponent ax ~ 1.3 for all the 
disorder configurations. 

Our experimental results can be compared now with 
the predictions of the nonlocal models [l7], |l8|, [n]] 
discussed in Section II. It appears that none of the mod- 
els can account for the different results obtained in the 
whole range of capillary numbers explored. Although the 
Flory-type argument used by Ganesan and Brenner Q 
is difficult to justify in this nonequilibrium situation, they 
obtain roughness exponents (ot\ — 3/4 and «2 = 1/2) 
consistent with our values for small and moderate Ca'. 
The exponents obtained by Hernandez-Machado et al. 
fl7f for the long length scales at moderate Ca', ol-i = 1/2 
and z = 1 (/?2 = 1/2), are in agreement with our ex- 
perimental values in the long time regime, where the 
viscous forces are dominant. The value of a for short 
length scales, corresponding to the short time regime of 
the model, gives a super-rough behaviour (ai > 1) domi- 
nated by surface tension in the plane, but larger than the 
exponent measured in the experiment. For large capil- 
lary numbers the numerical result of Dube et al. j^, ^| . 
ai = 1.3, is in agreement with our experimental results, 
but not the result /3i = 0.3 because, as pointed out above, 
these authors assume a spontaneous imbibition. At this 
point it would be interesting to know whether a model 
containing a quenched disorder both in the mobility and 
in the chemical potential would be able to explain the 
experimental results for the whole regime of Ca' stud- 
ied here, and specially the regime of large Ca' in forced 



imbibition. This remains an open question. 

Our last analysis is the variation of q c with v. The 
crossover wavenumber q c separates the regime in which 
the long length scale fluctuations (small q) are damped 
by the viscous pressure field, from the regime in which 
the short length scale fluctuations (large q) are damped 
by the interfacial tension in the plane of the cell. Since 
the relative importance of the viscous pressure field in- 
creases with Ca', it is expected that q c will increase with 
v. Specifically, a linear analysis of the interfacial prob- 
lem shows that the viscous damping is proportional to 
v\q\ and the interfacial tension damping is proportional 
to q 2 \q\, which results in q c oc u 5 jL7], [19] . 

In Fig. |22| we show the behaviour of q c with v for dif- 
ferent kinds of disorder. The error bars are relatively 
large because the exact location of q c in the experimen- 
tal power spectra is difficult to ascertain. For SQ 1.50 
we find that q c cx u ' 47 , in good agreement with the theo- 
retical prediction. It is interesting to note that, as we go 
to disorders of larger I (increasing persistence) , q c tends 
to be less sensitive to v. For T 1.50 we observe that q c 
becomes independent of v, within error bars. The rea- 
son can be understood in the framework of the following 
scenario: at a local level (see Fig. [j]) the motion of the 
interface in the SQ disorder can be viewed as a series 
of events formed by a period of nearly steady motion, 
a subsequent period of fast advance over a copper is- 
land (accompanied by an abrupt change of sign of the 
in-plane curvature), and a third period of fast advance 
over fiber-glass due to relaxation of the local in-plane 
curvature. Thus, damping of the short length scales due 
to the interfacial tension in the plane of the cell is basi- 
cally effective only in this last period, when the interface 
depins from the copper islands. As the persistence of the 
disorder in the y direction is larger, the relaxation peri- 
ods are less frequent. This can be observed through the 
histogram of Fig. [| When we increase the persistence 
of the disorder changing from SQ to SQ-n, the number 
of the smallest copper aggregations reduces almost one 
order of magnitude. In the limiting case of T 1.50 the 
disorder is continuous in the y direction, and the damp- 
ing role of the in-plane interfacial tension is effectively 
suppressed. 



VIII. CONCLUSIONS 

We have presented new experiments of forced fluid im- 
bibition in a Hole-Shaw cell with quenched disorder. We 
have used three main kinds of disorder patterns, SQ, 
SQ-n and T, characterized by an increasing persistence 
length I in the direction of growth. We have measured a 
robust roughness exponent [3 ~ 0.5 that is almost inde- 
pendent of the disorder configuration, interface velocity, 
and gap spacing, although the behaviour of the inter- 
facial width presents important fluctuations both dur- 
ing growth and at saturation, that progressively disap- 
pear as the disorder is more persistent in the y direc- 
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TABLE II: Summary of the values of the roughness exponents. 
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FIG. 22: Dependence of the crossover wavenumber q c on ve- 
locity v for three disorder configurations: SQ 1.50 (solid cir- 
cles), SQ-n 1.50 (open circles) and T 1.50 (triangles). 

tion. The roughness exponent a, however, shows a clear 
dependence on the experimental parameters, as summa- 



rized in Table y, discussed in the previous Section. Fi- 
nally, we have focused on the dependence of the crossover 
wavenumber as a function of the interface velocity v and 
I. For the shortest I, q c ~ v 0A7 , and becomes indepen- 
dent of v as the disorder is persistent in the y direction. 

The absence of pinning in experiments with disor- 
ders of increasing I allows a detailed investigation of the 
regime of large capillary forces. This regime presents 
some interesting novel features which are studied in de- 
tail in Ref. §§. 
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